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By using the Continuum RPA (CRPA) method, the incoherent transition strength of the exotic 
H~ — e~ conversion in the 208 Pb and 40 Ca nuclei is investigated. The question whether excited 
nuclear states lying high in the continuum give an important contribution to the incoherent rate 
is addressed. The admixture of spurious components in the rate coming from 1~ excitations is 
investigated in detail by using the self-consistent CRPA with Skyrme interactions as well as a less 
consistent version and by employing two ways to remove the spurious strength: the use of effective 
operators or simply the exclusion of the spurious state appearing close to zero energy. In all cases, 
the correction achieved is quite large. 



I. INTRODUCTION 

Theories that go beyond the Standard Model (SM) of 
particle physics allow for the violation of various sym- 
metries respected by the SM, among which the conser- 
vation of lepton flavour. Evidence for lepton flavour vi- 
olation (LVF) has been provided by the neutrino oscilla- 
tion experiments. In the charged-lepton sector, processes 
that would provide additional evidence for LVF and help 
distinguish among the various proposed mechanisms, in- 
clude the exotic neutrinoless conversion of a muon to an 
electron - where the neutrino and antineutrino involved 
in the conversion are assumed to be Majorana particles 
and can annihilate each other. 

In this context, the exotic conversion of a bound muon 
to an electron 
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has been studied both experimentally and theoreti- 
cally [JSSSI1SSSHQ3- The experimentally 
measured quantity is the ratio of the coherent rate, where 
the nucleus remains in its ground state, over the total 
capture rate. Although the coherent rate dominates the 
capture rate (accounting for about 90% of it), in order 
to calculate this ratio and make meaningful comparisons 
with experiment, both the coherent and the incoherent 
rate need to be evaluated theoretically. 

Various methods have been employed for the calcu- 
lation of the incoherent jjT — er conversion rate [J, 0, 
H 13 IH EH Elll 13 State-by-state calculations per- 
formed within the shell model and the quasiparticle RPA 
(QRPA) indicate that 

1. the main contribution to the incoherent rate comes 
from low-lying states 

2. the contribution of the 1~ channel is very large 
(about 50%, for all mechanisms leading to [i,~ — e~ 



conversion). Therefore, it is essential to properly 
remove possible spurious center-of-mass (CM) con- 
taminations. 

In a recent paper [l5j we used a Continuum-RPA 
(CRPA) method with Skyrme interactions to address the 
question as to how insignificant the contribution of ex- 
cited states lying high in the continuum really is. We 
considered natural-parity excitations of the 208 Pb nucleus 
and found that high-lying strength is not negligible. In 
this work we continue this investigation by looking at a 
lighter nucleus as well, namely Ca. We also examine in 
detail the admixture of spurious components in the rate 
coming from 1~ excitations. To this end we have used 
both the self-consistent CRPA and a non-consistent ver- 
sion (ncRPA) and have employed two ways to remove the 
spurious strength: by using effective dipole operators, as 
in Ref. |lq|. and by simply throwing away the spurious 
state appearing close to zero energy. In both cases, the 
correction achieved is quite large. 

Useful definitions and basic information on the meth- 
ods are provided in Sec. [HJ In Sec. IIHI we apply our 
CRPA method, which takes the full continuum into ac- 
count, to the nuclear targets 208 Pb and 40 Ca. The spu- 
rious strength is discussed in detail in Sec. IIIIBI We 
conclude in Sec. IIVI 



II. 



DEFINITIONS AND METHOD OF 
CALCULATION 



The inclusive (/i~,e~) rate is evaluated by summing 
the partial contribution of all final states \f). For spher- 
ical or nearly spherical nuclei, the vector contribution is 
given by |4| 



*■ = £ - 

f 



\(f\O a (q f )\0)\ 2 



(2) 



* Email panagiota. papakonstantinouOphysik. tu-darmstadt . de 



where O a (q/) represents the vector- type transition oper- 
ator resulting in the context of a given mechanism me- 
diated by a photon (a = 7), a JF-boson (a = W) or a 



Z-particle exchange (a — Z). Here g/, with magnitude 
qf = m^ — £(, — Ef, is the momentum transferred to the 
nucleus. Ef is the energy of the final state |/) with re- 
spect to the ground state |0), £& is the binding energy 
of the muon and m^ its mass. The transition operators 
have the form 



o a (q) = gvfv }2 6c "( r j) e iqr3 1 c «( t j) = I + i/ 3 ^ 

, (3) 

where Tj is the 3rd component of the jth particle's 
isospin. The parameter fy = 1.0 represents the vector 
static nucleon form factor and the normalization coef- 
ficient gv takes the value 1/6 for the photonic case and 
1/2 for the non-photonic W boson and SUSY Z exchange 
5j. The value of (3 a depends on the model assumed. We 
have adopted the values from Ref. @|. Thus, protons 
(neutrons) contribute to a given process with a "charge" 
whose value is determined by c a (l/2) = 1/2 + /3 ra /6 
(c a (— 1/2) = 1/2 — /3 a /6). In the photon and Z case, 
the isoscalar and isovector components of the transition 
operator are (almost) equally important, whereas Ow is 
predominantly isoscalar. 

By assuming that the initial and final states are of def- 
inite spin and parity, a multipole decomposition of the 
operators of Eq. J3J into operators T a ^M of orbital angu- 
lar momentum rank L can be carried out. For spherical 
nuclei we can assume, without loss of generality, q = z. 
Then, only terms with M = survive, for which we ob- 
tain {T aL (q) = T aLQ (q)) 



T aL (q) = gvfvVMZL + 1) J2 teaWJLiqrJYuofa). 

i=i 

(4) 
A phase factor (— i) L has been omitted. The contribution 
of each multipolarity to the transition rate 5* a reads 
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We now rewrite the rate S a L as the integral of a suitable 
distribution over excitation energy: 



S aL = JdER aL (E) 



with 



R aL (E) = 



E 1 
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In the above expression we have set k = 1 — tb/rn^, while 



The final states |/), excited by the single-particle op- 
erator T q l, are of particle-hole {ph) type. Then, the 
distribution R' aL (E), and from it R a L(E), can be calcu- 
lated following the standard RPA method. Subsequently, 
Eq. JBJ) can be used to evaluate the total rate S a L- 

We consider ph excitations, built on top of the mean- 
field ground state of a closed-shell nucleus and sub- 
jected to the ph residual interaction. In particular, the 
quantities introduced above are calculated using a self- 
consistent Skyrme-Hartree-Fock (SHF) plus Continuum- 
RPA (CRPA) model. The HF equations describing the 
ground state are derived variationally from the Skyrme 
energy functional. RPA excitations are considered on 
top of the HF ground state. The CRPA is formulated 
in coordinate space, so that the full particle continuum 
is taken into account. The same Skyrme interaction is 
used to calculate the ground state properties and the 
residual ph interaction. The model is described in detail 
in Ref. | 15| and refs. therein. In this work we have em- 
ployed the SkM*[l7| parametrization of the Skyrme force. 
It describes satisfactorily giant resonances of stable nu- 
clei, and therefore it is suitable for the present study. In 
order to test the sensitivity of our results on the inter- 
action used, we have also used MSk7 [13, which has a 
large effective mass, thereby shifting most excited states 
to lower energies compared to the more reliable SkM*. 

For the purposes of Sec. IIIIBI we have employed also 
the non-consistent CRPA version (ncRPA) of Ref. 19] 
and the corresponding numerical code |20j| . The ground 
state is described by a Woods-Saxon potential of radius 
(A — l) 1 / 3 ?^ (r = 1.25 fm) and diffuseness ao = 0.65 fm, 
including central (strength Vq = —53 MeV), spin-orbit 
(V so = 15.5 MeV fm 2 ), symmetry (V T = 20 MeV) 
and Coulomb terms. The residual ph interaction is 
a simplified Skyrme interaction without spin and ve- 
locity dependence (to = —1100 MeV fm 3 , xq = 0.5, 
£3 = 15000 MeV fm 6 ). Its strength is scaled by a fac- 
tor V^ ca i so as to to bring the spurious state close to zero 
energy. 

Results are presented and discussed in the next section. 



III. RESULTS 

Next, we present results for the nuclei 208 Pb and 40 Ca. 
The muon binding energy in 208 Pb is £b = 10.475 MeV. 
The particle threshold energy E t hr is 8.09 MeV in the 
case of the SkM* force. The respective values for 40 Ca 
are e b = 1.0533 MeV, E thr = 8.86 MeV. We have ob- 
tained results for L = 0, 1, . . . , 6 and for natural parity, 
(— 1) L . The most important contributions to the inco- 
herent transition rate are expected from L < 4 |lC|. 



R' aL (E) =J2\(f\ T ^(™»- £b- E f )\0)\ 2 6(E - E f ) (8) 
/ 

stands for the "strength distribution" corresponding to 
the operator T a L(q), with q = m^ — £& — E. 



A. Incoherent transition rate in the continuum 

In Fig. 2] first two panels, the distribution R a L(E) 
is plotted as a function of E, for the + and 4 + tran- 



sitions of 208 Pb. For the 1" distribution of 208 Pb the 
reader is referred to Fig. 0] (first panel, full lines) in 
the next subsection. In the monopole case, L = 0, the 
Isoscalar (IS) Giant Monopole Resonance (GMR) is the 
main peak. For 7 and Z, there is considerable contri- 
bution coming from higher energies (20-35 MeV), i.e., 
the isovector (IV) GMR region. For L = 1, the IV Gi- 
ant Dipole Resonance (GDR) corresponds to the strength 
clustered around E sa 12 MeV. For W and Z, important 
contribution seems to come from higher energies (above 
20 MeV), in particular, the IS GDR. For W exchange, the 
region below 10 MeV contributes significantly. In this re- 
gion we find the oscillation of the neutron skin against 
the nuclear core (pygmy dipole resonance) [l||. In the 
quadrupole case, L = 2 (not shown), both the IS Giant 
Quadrupole Resonance (GQR), close to 11 MeV, and the 
collective low-lying state are strong. There is some con- 
tribution from energies higher than 15 MeV, i.e., from 
the IV GQR region, especially in the cases 7 and Z. For 
L = 3 (not shown) the strength is mostly concentrated 
in the collective octupole state at low energy. For L > 3, 
eg. for L = 4, the calculated strength is quite fragmented 
and most of it lies below 20 MeV. 

In Fig. ^ last two panels, the distribution R a L{E) is 
plotted as a function of E, for the + and 3 _ transitions 
of 40 Ca. For the 1~~ distribution of 40 Ca the reader is 
referred to Fig. 0] (second panel, full lines) in the next 
subsection. The monopole strength distribution is dom- 
inated by the broad IS GMR. Above 27 MeV the IV 
GMR is excited through the 7 and Z mechanisms. The 
IV GDR dominates the dipole spectrum. In the 2 + distri- 
bution we find the strong IS GQR at about 17 MeV and 
seemingly little strength at higher energies. The 3~, 4 + 
(not shown) distributions are quite fragmented. A strong 
peak appears at about 5 MeV in the 5~ distribution (not 
shown) along with frargmented strength at higher ener- 
gies. In the 6 + distribution (not shown) we find a strong 
peak at around 15 MeV and a broad structure in the 
continuum between 20 and 30 MeV. 

In Fig. we plot the fraction of the total strength S a L 
coming from states below the particle threshold (S a L,thr) 
and the fraction coming from states below 20 MeV 
{SaL,20Mev), vs. the multipolarity L, for 208 Pb and using 
the Skyrme parametrizations SkM* and MSk7. We see 
that for low multipoles L = 0, 1,2 only a small portion 
of the strength originates from energies below particle 
threshold. The trend followed is similar for all mecha- 
nisms and independent of the interaction used. For even 
multipoles L — 2,4 a big portion of the contribution is 
pushed to higher energies as compared to the neighbor- 
ing odd ones. For some multipoles (L — for photonic 
mechanism, L = 1 for W^-boson exchange), a significant 
portion of the strength comes from above 20MeV. 

In Fig. |3 we show the fractions S a L,tht/S a L and 
S a L,20Mcv/S a L for both nuclei 208 Pb and 40 Ca, evalu- 
ated using the SkM* force. We also show, for 40 Ca, 
the fraction S a L,35Mev/SaL coming from states below 
35 MeV. The quantities S a L,thr/S a L for 40 Ca show that 



for L = 0, 1 and, in addition, for L = 2, 4, 6 almost all 
of the strength comes from above threshold. One should 
bare in mind that 40 Ca is an £— closed nucleus, with- 
out low-lying A7V — ph states. 208 Pb, on the other 
hand, is Ij— closed and low-lying L = 2, 4, 6 transi- 
tions between spin-orbit partners are present. For odd 
multipolarities in 40 Ca an important amount of strength 
comes from continuum excitations as well. A larger frac- 
tion of strength is found above 20 MeV for 40 Ca than 
for 208 Pb. In this sense, high-lying excitations are found 
more important for 40 Ca than for 208 Pb. Of course, the 
value of -E max = 20 MeV was chosen arbitrarily. The 
oscillator level spacing (hw ~ 41^4 _1 / 3 ) of a light nu- 
cleus is larger than for a heavier nucleus and the en- 
ergies of the corresponding excitations are higher. If we 
choose to compare, e.g., the S a L,2QMeV I S a L of 208 Pb with 
the S a L,35Mev/S a L of Ca (the maximum energy being 
« 3 shells in both cases), we find that, in the latter case, 
values closer to unity are reached. 

We have calculated also the fraction of the total 
strength S a L, coming from states below 50 MeV, for L 
up to 6 and for both nuclei. In all cases, the fraction 
is practically equal to unity. This means that the dis- 
cretized versions of RPA and QRPA are safe to use if the 
energy cutoff is large enough to sufficiently account for 
transitions below this value. 



B. Removal of spurious strength 

It is well known [13, uM that the 1~ excitations con- 
tain admixtures of the spurious excitation of the center of 
mass (CM) of the nucleus, corresponding to a situation in 
which the nucleus moves as a whole around the localized 
fictitious potential well. Normally, spurious components 
are separated out by the RPA methods. However, the use 
of a truncated model space and non-self-consistent single 
particle energies in ordinary RPA and the other versions 
of QRPA destroys the translational invariance and in- 
serts spurious excitations into the spectrum. Thus, the 
spurious CM state is not completely separated from the 
real (intrinsic) nuclear excitations, and in addition its en- 
ergy eigenvalue is not zero. In Continuum-RPA models 
with Skyrme interactions it has been possible to achieve 
a high degree of self-consistency, i.e. the same interaction 
is used for the HF calculation of ground state properties 
and for the residual interaction. In addition, no trunca- 
tion is involved. However, due to the formulation of the 
model in coordinate space, it is common practice to ex- 
clude the Coulomb and spin-orbit contribution (at least) 
to the residual interaction. Therefore, self-consistency is 
violated and, even in cases where the spurious state ap- 
pears very close to zero energy, some spurious strength 
may remain at higher energies. 

For electric dipole excitations, the problem is usually 
treated by using effective charges [22J. Similarly, in the 
case of IS dipole excitations, effective operators are used 
|23l 12 l| , which minimize the spurious admixture in the 
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The distribution R a L(E) in 208 Pb for L — 0, 4 and in 40 Ca for L — 0, 3. Skyrme parameterization SkM* has been 
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FIG. 2: Fraction of the total strength S a L, for the nucleus 
208 Pb, coming from states below the particle threshold (top) 
and below 20 MeV (bottom) vs the multipolarity L. Skyrme 
parameterizations SkM* and MSk7 have been used. Lines are 
drawn to guide the eye. 



FIG. 3: Fraction of the total strength S a L, for the nuclei 
208 Pb and 40 Ca, coming from states below the particle thresh- 
old (top) and below 20 MeV (bottom) vs the multipolarity L. 
For 40 Ca the fraction of strength coming from below 35 MeV 
is also shown (bottom). Skyrme parameterization SkM* has 
been used. Lines are drawn to guide the eye. 



strength distribution. (The effect on the IS dipole ex- 
citations of 208 Pb was examined in detail in Ref. |25|.) 
In Ref. |15J we presented a similar prescription for the 
operators involved in /i~ — e~ conversion. In particular, 
the operators T a \ which induce the 1~ excitations, 

A 
Tai(q) = 5>( r i)/fa) y i°ft);/( r ) = 6gvfvVl2^h(qr), 

(9) 
are replaced by respective effective operators 

A 

T£™(q) = EtcfaO/fa) - Var 3 ]Y w {f 3 ). (10) 

i=i 

In principle, the operators T al and T^°" induce the same 
intrinsic excitations, because they only differ by a term 
which translates the center of mass. The parameter r\ a 



in Eq. Q10[l is determined so as to eliminate the spurious 
CM excitation within the collective model. One finds [jji 



Va = gvfvlVSTTq 
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The point-proton and neutron form factors, F p (Q) and 
F n (Q) respectively, are calculated numerically using the 
ground-state densities. 

In Fig. 01 first two panels, we plot the dipole distri- 
butions R a i (E) of Eq. (jHJ , for photon, W- and Z-boson 
exchange diagrams, calculated by using the corrected and 
uncorrected operator. The SkM* interaction was used. 
Results are shown for both 208 Pb and 40 Ca. Most of 
the spurious strength below «6 MeV has been removed. 
The strength distributions above 20 MeV are practically 



unaffected. The strength in the region of the IVD reso- 
nance appeares redistributed. The effect of the correction 
appears strongest in the case of the W-boson exchange 
mechanism and for 40 Ca. The pygmy dipole state of 
208 Pb below 10 MeV is strongly affected by the correction 
in the photonic and W cases. 

In most of the cases, the effect of the correction above 
the tail of the spurious state appears small. There- 
fore, it would have been a fairly good approximation to 
simply throw away the spurious state and calculate the 
transition rate coming from states above, approximately, 
6 MeV. Such an approximation was used in Ref. J10| . 
One should keep in mind, however, that our CRPA model 
is self-consistent to a high degree. Thus, the excitation 
spectrum above the spurious state is almost free of spuri- 
ous components, already before the operator correction, 
and therefore almost insensitive (up to 10 %, as we will 
see) to the change of operators. This is not necessar- 
ily the case in less consistent models such as the one in 
Ref. [il| and the ncRPA model described in Sec. [H] 

Next we employ the ncRPA method and calculate 
again the R a \(E) distributions using the operators T a \ 
and !X r - In the case of 208 Pb ( 40 Ca) the residual inter- 
action was scaled by a factor V sca \ = 0.94 (1.004) in order 
to bring the spurious state close to zero energy. The re- 
sults are shown in Fig. 0J last two panels. Again, in all 
cases, most of the strength of the spurious state is elim- 
inated. In general, the effect appears small for 208 Pb, 
above the tail of the spurious state. The pygmy dipole 
state is strongly affected and the strength of the IS GDR 
is slightly reduced (see W— boson mechanism). For 40 Ca 
the effect appears small in the region of the IV GDR, but 
not below and above. In particular, the use of corrected 
operators results in significant additional strength above 
20 MeV. In total, the transition strength above the spuri- 
ous state (above 6 MeV) increases. The sensitivity of the 
calculation on the operators used implies a bad degree of 
consistency which renders the calculation unreliable. 

In Tabled we list, for both nuclei, i) the portion s s ^ t 
of transition strength removed from the total contami- 
nated 1~ transition strength S a i when corrected opera- 
tors T£°™ are employed, ii) the portion of strength S sp 
carried by the spurious state when the uncorrected T a \ 
operators are used (calculated as the portion of strength 
lying below 6 MeV) and iii) the portion s!f 6MGV of the 
strength removed from above 6 MeV excitation energy 
when corrected operators are used (with respect to the 
uncorrected strength above 6 MeV). Results have been 
obtained with both models, i.e., our self-consistent CRPA 
model with interaction SkM* and the non-self consistent 
CRPA model labelled ncRPA. 

Let us first examine the behavior of the self-consistent 
CRPA model. As the values of s s t ^ t show, about 90% of 
the total transition rate was spurious in all cases. We ex- 
pect this result to be independent of the nucleus and the 
interaction used, because the spurious state at low energy 
always dominates the isoscalar dipole strength distribu- 
tion (which contributes in all three mechanisms) and be- 



cause the corrected operators are, by construction, most 
effective for this state (thus removing practically all its 
strength). We were not able to demonstrate this in the 
particular case of 40 Ca, because the energy of the spu- 
rious state in this case was found imaginary. In other 
words, we were not able to evaluate and take into ac- 
count properly the strength of the spurious state, before 
or after the correction. Therefore, the respective num- 
bers in Tableware placed in brackets. 

Suppose now that, in order to evaluate the rate coming 
from intrinsic 1~ excitations, we would apply the proce- 
dure of using the T a \ operators and then simply excluding 
the strength of the spurious state from calculating the to- 
tal rate. As the values of S sp show, for 208 Pb, we would 
have removed roughly as large a portion as s t „ t . Again, 
for 40 Ca it is not possible to conclude. 

From Table [3 we notice that for stf 6Mc y and for a = 7 
the result is small in absolute value, but negative, for 
208 Pb. It represents the numerical accuracy of our cal- 
culation and, being small, it indicates that the degree of 
self consistency reached by our HF+CRPA model is suf- 
ficient to achieve a satisfactory separation of the spurious 
transition in this case. For the W and Z cases, however, 
spurious admixtures of more than 6% are found above 
6 MeV. These numbers (which are not free of numer- 
ical inaccuracies, as explained in Ref. [l5|]) vary when 
different Skyrme interactions are used, with their values 
remaining below 10%. In the case of 40 Ca, the values of 
s >6MoV f° r ^ ne T anc ^ ^ mechanisms are lower than 10% 
as well. Its value for W-exchange is much larger. Most of 
it, however, comes from below 15 MeV (below the GDR 
region) . 

When the ncRPA model is used, similar results are ob- 
tained for s^Q t and S sp as with the self-consistent model. 
The 208 Pb excitation spectrum above the tail of the spu- 
rious state is a little more contaminated than within the 
self-consistent CRPA, as the s^ 6Mc v values show. As for 
the 40 Ca nucleus, when corrected operators were used, 
the strength above 6 MeV was significantly increased. 
Thus, the corresponding values of s!f 6McV are large and 
negative, indicating a bad performance of the model. 

From our results, obtained for two nuclei within two 
different continuum- RPA models, one may conclude that 
more than 85% of the total 1~ rate, when no correction is 
considered, is spurious. Whether effective operators are 
used to remove the total spurious strength, or whether 
one excludes the spurious state from the calculation of 
the total 1 _ , the above conclusion remains the same. 
This does not mean that there are no spurious contam- 
inations above the zero-energy spurious mode. When 
making use of corrected operators, the strength of the 
excited states across the spectrum is redistributed. The 
overall effect, however, is to reduce the total S a i rate 
by the amount that was initially carried by the spurious 
state. 

Within the QRPA calculations of Ref. [l(| , for six nu- 
clei (including 208 Pb), a more moderate correction (less 
than 60%) on S a i was achieved by considering the lowest 
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FIG. 4: The dipole distributions R a i(E), for 7-photon and W-boson exchange diagrams of the fi~ — > e~ conversion in 208 Pb 
and 40 Ca. The results have been calculated for the dipole operator T a i (dotted line) and for the corresponding corrected 
operator given by Eqs. ill 01 . pip . First two panels: The self-consistent CRPA with Skyrme parametrization SkM* was used. 
Last two panels: the non-self-consistent RPA was used. 



1 _ state as purely spurious and simply excluding it. One 
should keep in mind the differences between the models 
used there (QRPA calculations with a renormalized G- 
matrix interaction) and the ones used here. First of all, 
priority was given in Ref. [l(j to the best possible repro- 
duction of experimental spectra, rather than controlling 
the self-consistency. In addition, a truncated model space 
was used, contrary to the methods used in the present 
work. The various approximations entering may have in- 
fluenced the spuriosity results in an unpredictable way. 
Both consistency and completeness of the space are im- 
portant in order to move as much spurious strength as 
possible close to zero energy. Of course, the lowest V 
state appeared very close to zero energy. It was found 
that it was spurious by 60-80% (except for 126 Yb, for 
which it was about 90% spurious). Thus, the approxi- 
mation of considering it as purely spurious is not a very 
good one. The numbers mentioned above refer to the 
overlap of the lowest 1 _ state with the "purely spurious" 
RPA state. Approximations were inevitable when con- 
structing the purely spurious state, in order to normalize 
it. 



IV. CONCLUSIONS 
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9.1 




40 Ca- 


ncRPA 


incr. 


incr. 


incr. 



TABLE I: Percentage of the total 1 transition strength S a \ 



(s^ t ) and of the strength above 6 MeV (s : 



sp 

>6McV 



) consumed by 



spurious transitions, and percentage of strength below 6 MeV 
(~ strength of the spurious state, S ap ) for the three channels 
7, W, Z. Results are shown for the nuclei 208 Pb and 40 Ca 
obtained by using two models: the self-consistent Continuum 
RPA with Skyrme force SkM* and the non-consistent RPA 
(ncRPA). Last row: When corrected operators are used in 
this case, the strength above 6 MeV increases significantly. 
Square brackets: the energy of the spurious state is found 
imaginary, therefore its strength is not properly evaluated. 



We have investigated the incoherent rate of the exotic 
H~ — e~ conversion in the nuclei 208 Pb and 40 Ca. We em- 
ployed the Continuum- RPA method which is appropriate 
for explicit construction of the excited states lying in the 
continuum spectrum of the nuclear target. We used a 
self-consistent CRPA and Skyrme interactions to investi- 
gate the transition strength coming from natural-parity 
ph excitations up to L = 4. We found that a significant 
portion of the incoherent \x~ — e~ rate comes from high- 
lying nuclear excitations. The admixture of spurious 
components in the rate coming from 1~ excitations was 
investigated in detail by using the self-consistent CRPA 
with Skyrme interactions as well as a less consistent ver- 



sion of CRPA and by employing two ways to remove the 
spurious strength: the use of effective operators or simply 
the exclusion of the spurious state appearing close to zero 
energy. In all cases, we found that the greatest portion 
of the 1~ transition strength is due to the spurious CM 
excitation, a result in agreement with that of an exact 
method constructed recently yjj for removing spurious 
cont aminat ions . 
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